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Abstract: We investigate the strong-coupling phases that may arise in 3D Dirac and Weyl 
semimetals under the effect of the long-range Coulomb interaction, considering the many- 
body theory of these electron systems as a variant of the conventional fully relativistic 
Quantum Electrodynamics (QED). Eor this purpose, we apply two different nonperturba- 
tive approaches, consisting in the sum of ladder diagrams and taking the limit of a large 
number N of fermion flavors. We benefit from the renormalizability that the theory shows 
in both cases to compute the anomalous scaling dimensions of different operators exclu¬ 
sively in terms of the renormalized coupling constant, allowing us to determine the precise 
location of the singularities signaling the onset of the strong-coupling phases. We show 
then that the QED of 3D Dirac semimetals has two competing effects at strong coupling. 
One of them is the tendency to chiral symmetry breaking and dynamical mass generation, 
which are analogous to the same phenomena arising in the conventional QED at strong 
coupling. This trend is however outweighed by the strong suppression of electron quasipar¬ 
ticles that takes place at large N, leading to a different type of critical point at sufficiently 
large interaction strength, shared also by the 3D Weyl semimetals. Overall, the phase dia¬ 
gram of the 3D Dirac semimetals turns out to be richer than that of their 2D counterparts, 
displaying a transition to a phase with non-Eermi liquid behavior which may be observed 
in materials hosting a sufficiently large number of Dirac or Weyl points. 
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1 Introduction 


In recent years, condensed matter physics has witnessed several remarkable discoveries that, 
starting from graphene[l], have shown the existence of electron systems where the quasipar¬ 
ticles have linear dispersion about degenerate points at the intersection between valence 
and conduction bands. This has led to introduce the class of so-called Dirac semimet¬ 
als, in which the low-energy excitations can be described in terms of a number of Dirac 
spinors. Thus, there is already clear evidence that materials like NasBi or Cd 3 As 2 have 
Dirac points in their electronic dispersion [2-4], providing in that respect a 3D electronic 
analogue of graphene. More recently, there have been claims that materials like TaAs or 
the pyrochlore iridates should be examples of 3D Weyl semimetals [5, 6], characterized by 
having Weyl points hosting each of them fermions with a given chirality. 

Since the appearance of graphene, the possibility that the quasiparticles of a condensed 
matter system may behave as Dirac fermions has been very appealing, as the particular 
geometrical features of these mathematical objects have shown to be indeed at the origin of 
several unconventional effects. The Klein paradox[7] and the difficulty that such quasiparti¬ 
cles have to undergo backscattering[8] are a direct consequence of the additional pseudospin 
degree of freedom inherent to Dirac fermions. The topological insulators and their peculiar 
surface states can be taken also as an illustration of the fruitful interplay arising between 
the symmetry properties of Dirac fermions and the topology of the space[9, 10]. 

Another important feature of known Dirac semimetals is that they need to be modeled 
as interacting electron systems which are naturally placed in the strong-coupling regime. 
This is so as the long-range Coulomb repulsion becomes the dominant interaction, while 
the Fermi velocity of the electrons takes values which are at least two orders of magnitude 
below the speed of light. Thus, the equivalent of the fine structure constant in the electron 
system can be more than one hundred times larger than its counterpart in the conventional 
fully relativistic Quantum Electrodynamics (QED). This leads to expect relevant many- 
body effects as, for instance, the imperfect screening of impurities with sufficiently large 
charge on top of graphene[ll-14]. Another important property is the scale dependence of 
physical observables, like that predicted theoretically for the Fermi velocity in 2D Dirac 
semimetals [15, 16] and measured recently in suspended graphene samples at very low dop¬ 
ing levels[17]. We may think therefore of the Dirac semimetals as an ideal playground to 
observe scaling phenomena that would be otherwise confined to the investigation of field 
theories in particle physics. 

In any event, a proper account of the scaling effects must rely on the renormalizability 
of the many-body theory. The formulation of the interacting theory of Dirac semimetals 
requires the introduction of a high-energy cutoff in the electronic spectrum, which can 
be set to a rather arbitrary level. In order to ensure the predictability of the theory, 
it becomes therefore crucial that all the dependences on the high-energy cutoff can be 
absorbed into the redefinition of a finite number of parameters. In the case of Dirac 
semimetals, such a renormalizability cannot be taken for granted, as the many-body theory 
does not enjoy the full covariance which enforces that property in typical relativistic field 
theories. Nevertheless, the investigation of the 2D Dirac semimetals has already provided 
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evidence that the interacting field theory in these electron systems is renormalizable[18]. 
This condition has allowed for instance to compute very high order contributions to the 
renormalization of relevant order parameters, leading to an estimate of the critical coupling 
for chiral symmetry breaking that agrees with the value obtained with a quite different 
nonperturbative method like the resolution of the gap equation[19]. 

The aim of this paper is to investigate the strong-coupling phases that may arise in 
3D Dirac and Weyl semimetals under the effect of the long-range Coulomb interaction. In 
this respect, the many-body theory of these electron systems can be viewed as a variant of 
conventional QED. As already pointed out, a main difference with this theory lies in that 
the Fermi velocity of the electrons does not coincide with the speed of light, which makes 
that quasiparticle parameter susceptible of being renormalized and therefore dependent 
on the energy scale. In three spatial dimensions, the electron charge needs also to be 
redefined to account for cutoff dependences of the bare theory, which renders the question 
of the renormalizability even more interesting in the present context. 

In this paper we will apply two different nonperturbative approximations to the QED of 
3D Dirac and Weyl semimetals, namely the ladder approach and the large-approximation 
{N being the number of different fermion flavors), which can be viewed as complementary 
computational methods. In both cases, we will see that the field theory appears to be 
renormalizable, making possible to absorb all the dependences on the high-energy cutoff 
into a finite number of renormalization factors given only in terms of the renormalized 
coupling. 

We will show that the QED of 3D Dirac semimetals has two competing effects at strong 
coupling. One of them is the tendency to chiral symmetry breaking and dynamical mass 
generation, which are analogous to the same phenomena arising in the conventional QED 
at strong coupling[20-27]. This trend is however outweighed by the strong suppression of 
electron quasiparticles that takes place at large N, leading then to a different type of critical 
point at sufficiently large interaction strength[28], shared also by the 3D Weyl semimetals. 
We will see that the nonperturbative approaches applied in the paper afford a precise 
characterization of the respective critical behaviors in terms of the anomalous dimensions 
of relevant operators. At the end, the phase diagram of the 3D Dirac semimetals turns 
out to be richer than that of their 2D counterparts, displaying a transition to a phase 
with non-Fermi liquid behavior (genuine also of the 3D Weyl semimetals) which may be 
observed in materials hosting a sufficiently large number of Dirac or Weyl points. 

2 Field theory of 3D Dirac and Weyl semimetals 

Our starting point is the field theory describing the 3D fermions with linear dependence 
of the energy £ on momentum p, e(p) = u_f|p|; and interacting through the scalar part 
(p of the electromagnetic potential. The fermion modes can be encoded into a number N 
of spinor fields "ipi, i = 1,... where the index i represents in general the way in which 
the modes are distributed into different Dirac or Weyl points in momentum space. For 
convenience, we will choose to pair the fermion chiralities into four-component spinors, 
extending the notation to think of each ijji field as being made of two different chiralities 
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from a given Dirac point, or from two different Weyl points in the case of a Weyl semimetal. 
The hamiltonian can be written then as 


/ 


H = -ivF / d^r V'i(r)7 ' + eo / d^r V'i(r) 7 oV’i(r) 4>{r) 


/' 


( 2 . 1 ) 


where = il’j'jo and { 70 -} is a collection of four-dimensional matrices such that { 7 ^, 7 ,^} = 
2 diag(l,-l,-l,-l). 

In order to complete the field theory, we have to provide the propagator of the scalar 
field (j), for which we take the Lorentz gauge in the dynamics of the full electromagnetic 
potential. Thus, we get from the original relativistic theory 





^ iq-(r-r') 
27r 
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ijj^ jc? — ir] 
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We are interested in systems where the Fermi velocity Ui? is much smaller than the speed of 
light c, which makes possible to adopt a simpler description taking the limit c ^ 00 . This 
justifies that we can neglect magnetic interactions between our nonrelativistic fermions, 
and that we can just deal with the zero-frequency limit of (2.2). Thus, the free propagator 
for the scalar potential will be taken in what follows as 

^o(q) = 4 (2-3) 


A remarkable property is that the hamiltonian (2.1) together with the interaction 
mediated by (2.3) define an interacting field theory that is scale invariant at the classical 
level. That is, under a change in the scale of the space and time variables 


t —)• st , r —)■ sr 


(2.4) 


the helds can be taken to transform accordingly as 

ct){r) ■^())(r) , Tpir) ^V’(r) (2.5) 

This gives rise to a homogeneous scaling of the hamiltonian (2.1), consistent with the 
transformation of an energy variable. 

Such a scale invariance has important consequences, since it makes possible to give 
meaning to the divergences that the held theory develops at high energies. The compu¬ 
tation of the quantum corrections requires indeed the introduction of a high-energy cutoff 
that spoils the classical scale invariance. However, provided the held theory is renormaliz- 
able, the singular dependences on the cutoff can be still absorbed into a hnite number of 
renormalization factors redehning the parameters of the theory. Writing the action S cor¬ 
responding to the hamiltonian ( 2 . 1 ), we may start with a formulation introducing suitable 
renormalization factors Z^,Z^,Ze and 

S = j dtd^r Z^'il;i{r){i'yodt + iZyVR-f-d)'il;i{r)-Zee J dtd^r Z^ipi{r)'yo'ilji{r) Z^(j){r) ( 2 . 6 ) 
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The gauge invariance of the model implies that Zf.Z^ = 1. The point is that, if the field 
theory is well-behaved, one has to be able to render all the physical observables cutoff- 
independent, when written in terms of the renormalized parameters vr and e. 

The first example of the need to implement a high-energy cutoff can be taken from the 
computation of the electron self-energy. The first perturbative contribution is given by the 
first rainbow diagram at the right-hand-side in Fig. 1, which corresponds to the expression 


zi:i(k) = -el 



du}p -Joujp + VFl • p 1 

ivr — IT] (k — pY 


(2.7) 


This contribution to the electron self-energy does not depend on frequency, but shows a 
divergent behavior at the upper end of the momentum integration, pointing at the renor¬ 
malization of the Fermi velocity vp. In what follows, instead of dealing with a cutoff in 
momentum space, we will prefer to use a regularization method that is able to preserve the 
gauge invariance of the model. For this purpose, we will adopt the analytic continuation 
of the momentum integrals to dimension D = 3 — e, which will convert all the high-energy 
divergences into poles in the e parameter [29]. 

In the dimensional regularization method, the bare charge eg must get dimensions 
through an auxiliary momentum scale fi, being related to the physical dimensionless charge 
e by 

eo = (2.8) 

The computation of the first-order contribution to the self-energy (2.7) gives then 


Si(k) 


d^p 1 1 

2^7 (27r)'^"^ *^|p| (k - p)2 

^ f d^p (1 — 

4^7 ^ Jq *[(k-p)2a;-Fp2(l- x)]3/2 

^ J (2^ Jo "^'^[p2+k2x(l-x)]3/2 

["dx 

20F^ (47r)^/2 ^ Jo [k2x(l - x)]3/2-^/2 

F(i6)F(2-|)r(l-f) 
( 47 r) 2 -e/ 2 T |]^|, r(|-e) 


(2.9) 


The expression (2.9) shows the development of a 1/e pole as e —)• 0. This divergence 
can be indeed absorbed into a renormalization of the Fermi velocity as the full fermion 
propagator G(k, w) is related to the self-energy E(k,a;) by 

G(k,a;)“^ = ^^(yow - ZyVRj ■ k) - Z^S(k,a;) (2.10) 


Thus, to first order 
taking 


in e^, the fermion propagator can be made finite in the limit e —0 by 


e 


2 
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Zy — 1 — 


Gtt'^vr e 


( 2 . 11 ) 

















This renormalization of the fermion propagator already gives an idea of the effective 
dependence of the Fermi velocity on the energy scale. The bare Fermi velocity Z^vr 
cannot depend on the momentum scale since the original theory does not know about 
that auxiliary scale. We have therefore 

^^{Z,vr) = ^ (2.12) 


For the same reason, the bare coupling eo cannot depend on which leads to the equation 


d 






(2.13) 


Combining (2.12) and (2.13), we arrive at the scaling equation 


d 1 

on OTT^ 


2 


(2.14) 


This is the expression of the growth of the renormalized Fermi velocity vr in the low- 
energy limit, approached here as /i —)• 0. This scaling, which parallels the behavior of the 
Fermi velocity in 2D Dirac semimetals, is the main physical property deriving from the 
lowest-order renormalization of the 3D Dirac and Weyl semimetals [30]. 

To close this section, we comment that the result (2.9) for the electron self-energy has 
actually a wider range of validity, beyond the lowest-order approximation in which it has 
been obtained. This can be seen by exploiting a feature that is also shared with the 2D Dirac 
semimetals in the so-called ladder approximation, that is when corrections to the interaction 
vertex and the interaction propagator are neglected in the electron self-energy [18]. Then 
we remain with the sum of diagrams encoded in the self-consistent equation represented in 
Fig. 1. The electron self-energy in this approximation, Siadder(k)) niust have the form 


Sladder(k) = /(k) 7 • k 


(2.15) 


and therefore it is bound to satisfy the equation 


^^ladder(k) 


iSi(k) eo 


f d^p f dup^ a;2 + ^;2,p2 

J (27r)^y 2vr^'"''"'^''^P^(-a;2 + 4p2-*^)Mk-p)2 


(2.16) 


It is now easy to see that the second term at the right-hand-side of (2.16) identically 
vanishes. By performing a Wick rotation Up = ioJp, we get 

r dujp ujI + . r dujp -^1 + ^ .2 171 

J 2 tt {- 0 }“^ + - ir])^ V 27r (aj2+^2^p2)2 

showing that 

5^1adder(k) = Fli(k) (2-18) 

The result (2.18) has important consequences for the ladder approximation, since it 
implies that the low-energy scaling of the Fermi velocity is again the main physical effect 
derived from the electron self-energy, limiting the possible electronic instabilities that can 
arise in that nonperturbative approach. 
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Figure 1. Diagrammatic representation of the ladder approximation for the electron self-energy. 

3 3D Dirac semimetals in ladder approximation 

The iterative sequence of interactions, illustrated in Fig. 1 for the electron self-energy, 
defines a partial sum of perturbation theory which can be also applied to analyze the 
renormalization of composite operators. This provides a suitable approach to investigate 
possible instabilities of the electron system, since some of those operators correspond to 
order parameters characterizing the breakdown of respective symmetries of the field theory. 
Once the renormalization program is accomplished, the singularities found in the correlators 
of the composite operators may signal the condensation of a given order parameter, pointing 
at the onset of a new phase of the electron system. 

The most relevant composite operators are made of bilinears of the fermion fields, 
having associated vertex functions that can be used to measure different susceptibilities of 
the electron system. At this point, the discussion is confined to the 3D Dirac semimetals 
(excluding the case of Weyl semimetals), for which we can write different composite oper¬ 
ators in terms of four-component spinors about the same Dirac point in momentum space. 
We pay attention in particular to those operators corresponding to the charge, the current. 


and the fermion mass, given by 

Po(r) = V'iWToV'iW (3.1) 

Pc(r) = V’i(r)7^i(r) (3.2) 

Pm{r) = '0i(r)V'i(r) (3.3) 

The respective one-particle-irreducible (IPI) vertices are 

ro(q,Wg;k,a;fc) = (po(q, Wq)V'i(k-F q, w*,-h Wq)V'i(k, Wfc))ipi (3.4) 

rc(q,a;g;k,a;fc) = (Pc(q> + q, Wfc + Wg)'0i(k, Wfc))ipi (3.5) 

= (/9m(q,^^g)V’i(k + q,Wfc+ a;g)V>i(k,Wfe))iPi (3.6) 


The point is that, assuming the renormalizability of the field theory, there must exist 
a multiplicative renormalization of the vertices 

TOjren — -^0^0 
Tc,ren — 

Prrijren — 

that renders To^renj Tc,ren and Tm.ren independent of the high-energy cutoff. We check 
next that property of the field theory, dealing with the same iteration of the interaction 


(3.7) 

(3.8) 

(3.9) 
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Figure 2. Self-consistent diagrammatic equation for a generic vertex in the ladder approxima¬ 
tion. 

applied before to the electron self-energy. In the present context, this defines the ladder 
approximation as given by the self-consistent diagrammatic equation represented for a 
generic vertex in Fig. 2. 

We first observe that the renormalization of Fq must be dictated by the gauge invari¬ 
ance of the model, since the composite operator /9o(r) already appears in the action ( 2 . 6 ). 
This implies the result 

^0 = (3.10) 

meaning that, according to the previous analysis of the self-energy, it must he Zq = 1 
in the ladder approximation. Moreover, the composite operator may be introduced in 
the action multiplying it by the vector potential, showing that its renormalization can be 
related by gauge invariance to that of the kinetic term in (2.1). We conclude therefore that 

Zc = Z^Z^ (3-11) 

which, in the ladder approximation, leads to Zc = Z^. 

The results (3.10) and (3.11) can be also obtained more formally from the Ward iden- 
dities that reflect the gauge invariance of the model at the quantum level. They lead in 
particular to the equations 

^ ( 70 a; - S(k, uj)) = ro(0,0; k, ui) (3.12) 

• k-h S(k,a;)) = rc(0,0;k,a;) (3.13) 

Vp OK 

These identities were used in Ref. [18] to check the suitability of the dimensional regular¬ 
ization method to preserve the gauge invariance of the field theory of 2D Dirac semimetals. 
It was shown there that, if the the electron self-energy is computed with the set of diagrams 
encoded in the equation of Fig. 1, then the Ward identities are satisfied when one takes 
for the vertex the series of ladder diagrams, but dressed with the lowest-order electron 
self-energy correction. With this scope of the nonperturbative approach, it can be also 
shown by direct renormalization of Tc that Zc becomes as simple as the expression given 
by ( 2 . 11 ), which arises then from a remarkable cancellation of different contributions in 
the computation of the vertex. 

We conclude that the charge operator is not renormalized, while the current operator 
requires a multiplicative redefinition by the simple factor ( 2 . 11 ) in the ladder approxima¬ 
tion. The absence of a singularity for any value of the coupling constant excludes the 
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possibility of having an instability related to the condensation of the charge or the cur¬ 
rent in the electron system. We turn now to the case of the mass vertex (3.6), whose 
renormalization is not dictated by the previous analysis of the electron self-energy. 


3.1 Mass vertex and dynamical mass generation 

Our starting point for the analysis of the mass vertex in the ladder approximation is the 
self-consistent equation represented in Fig. 2. We are going to be mainly interested in 
the renormalization of F^ and, for that purpose, it is enough to consider the limit of 
momentum transfer q —)• 0 and Uq —>■ 0. The vertex must satisfy then the equation 


Fm( 0 , 0 ;k,a;fc) = 

d^p (hjp -7oWp -I- vpl ■ p 


1 + iCg 


(27r)'^ 27r —a;| -|- — irj 


^ -70Wp -h VFl ' P 

—Fm(0,0; p, ^p)—:r^j2Z3. -“70- 


1 


—ojji + u|,p2 — irf (k — p) 


(3.14) 


''p I ^'1 

The resolution of (3.14) implies that F^ has to be proportional to the identity ma¬ 
trix. Moreover, it turns out to be a function independent of the frequency in this 
approximation. After a little algebra, we arrive then at the simplified equation 


rm(0,0; k, Wfe) ~ ^ ^ ~2 


d^p 


Fm( 0 , 0 ;p,Wfc) 


1 


1 


(3.15) 


/o \n-^ I I /I \9 

Eq. (3.15) can be solved by means of an iterative procedure, expressing the vertex as 
a power series in the effective coupling Ag = eg/dvrui?. 


F™(0,0;k,a;A,) = l + ^A^ FW(k) 


(3.16) 


n=l 


Assuming that the term Fm^(k) is proportional to l/|k|"''', the next order can be computed 
consistently from (3.15), taking into account that 


p(n-l-l) 


(k) ~ 1 / 


d^p 


(27r)-^ |pP+’^^ (k — p)2 


(dvr 




Ihvr^/^ p ^ l+ne "^ p ^ 3—{n+2)€ 

Thus, the vertex can be written as an expansion 

CXD 

F™(0,0;k,u;fc) = l +J]A^^ 


|k|(’T-+l)e 


(3.17) 


(3.18) 


n=l 


where each order can be obtained from the previous one according to the recurrence relation 


with 


■Sn+l — -dn+l(c) Sn 




( 4 ,r)-/^r(!aie)r(i- 
40r p(i^)p 


(n^) p (Ip) 

^ 3-(n-H2)£ ^ 


(3.19) 


(3.20) 
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We observe from (3.18) and (3.19) that develops higher order poles in the e param¬ 
eter as one progresses in the perturbative expansion. At this point, one has to check the 
renormalizability of the theory by ensuring that all the poles can be reabsorbed by means 
of a redefinition of the vertex like that in (3.9). In terms of the dimensionless coupling 

dvrui? 

the renormalization factor Zm must have the structure 

= (3.22) 

i=l 


(3.21) 


with residues di depending only on A. 

In the present case, it may be actually seen that the renormalized vertex rm,ren can 
be made finite in the limit e —)• 0, with an appropriate choice of the functions dj(A). We 
obtain for instance for the first perturbative orders 


di(A) 


d2{X) 

dsW 

d^X) 

d5{X) 

deW 


TT^ 


TT 


,3 36. 


IOOtt^ 


+ 7r'^) a' 


; 1 


- TT^) A" 


TT® 


/ Tvr^ 23vrn e 
V" 9 + 1440j ^ +••• 


^ (57 - 2vr2) A^ + ^ (lI4 - Tvr^) A^ 
6340 - 20I007r2 ISTvr'^) A® -b ... 

-^-^{33-7:^) 

4 727^5 I vr6 


_ vr® \ 48 

' + (75 - 27r2) A® - 

4327r6 ^ 

r A® + ... 


18, 


A® + .. 




(3.23) 

(3.24) 

(3.25) 

(3.26) 

(3.27) 

(3.28) 


An important point about the residues di{X) is that they can be chosen without having any 
dependence on the momentum k of the vertex. This is a signature of the renormalizability 
of the theory, by which we are able to render it independent of the high-energy cutoff, 
redefining just a finite number of local operators. 

On the other hand, the renormalization factor contains important information 
about the behavior of the theory in the low-energy limit. This stems from the anomalous 
scale dependence that the vertex gets as a consequence of its renormahzation[31]. The 
bare unrenormalized theory does not know about the momentum scale /r, but the factor 
Zjn lends to rm,ren Sin anomalous scaling of the form 




7m 


(3.29) 


- 10 - 

















(3.30) 


The anomalous dimension 7 ^ can be thus obtained as 

A* dZjYi 


The renormalization factor Zm is given by an infinite series of poles in the e parameter, 
and then it is highly nontrivial that the computation of 7 ^, from (3.30) may provide a finite 
result in the limit e —)• 0. At this stage, the dependence of Zm on /x arises from the scaling 
(2.13), since no self-energy corrections are taken into account yet. That leads to the 
equation 




<iiW 


00 , . 


7m — ^2-^ (=* (3.31) 

i=l / i=0 

The term with no poles in the e parameter already gives the result for the anomalous 
dimension 

d 

(3.32) 


7m = -A—di(A) 


But it remains however to ensure that the rest of the poles identically cancel in Eq. (3.31), 
which requires the fulfillment of the hierarchy of equations 


^di+i(A) = di(A)^di(A) (3.33) 

dX dX 

Quite remarkably, we have checked that the set of equations (3.33) is satished, up to 
the order we have been able to compute the residues di{X) in the ladder approximation. In 
this task, we have managed to obtain numerically the first residues up to order A^®, with a 
precision of 36 digits. Besides verifying the hierarchy (3.33) at this level of approximation, 
we have also analyzed the trend of the power series for such large perturbative orders, find¬ 
ing that a singularity must exist at a certain critical coupling. Concentrating in particular 
on the function di(A), we have found that the terms in its perturbative expansion 


di(A) = ^dS”U” (3.34) 

n=l 


approach a geometric sequence at large n. The values we have obtained for the coefficients 

in) 

d\ are represented in Fig. 3 up to order n = 28. Prom the precise numerical computation 
of the coefficients, it is actually possible to verify that their ratio follows very accurately 
the asymptotic behavior 


dl") , d’ d" d'" 

-7T — d H-1- 7T H-T + 


(3.35) 


From the estimate of the limit value d, we have obtained the finite radius of convergence 
for a value of the coupling constant 


Ac = ^ « 1.27324 (3.36) 

The coupling Ac corresponds to a point where the anomalous dimension 7 ^ diverges, 
which signals in turn the development of a singular susceptibility with respect to the mass 
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Figure 3. Plot of the absolute value of the coefhcients in the expansion of di(A) as a power 
series of the renormalized coupling A. 


parameter. That is, Ac has to be viewed as a critical coupling above which the electron 
system enters a new phase with dynamical generation of mass [32]. This characterization 
is also supported by the fact that, in the case of the 2D Dirac semimetals, a similar renor¬ 
malization in the ladder approximation [18] has shown to lead to a value of the critical 
coupling that coincides very precisely with the point for dynamical mass generation de¬ 
termined from the resolution of the mass gap equation[19]. It is indeed remarkable that 
these quite different approaches give an identical result for the point of chiral symmetry 
breaking in the electron system. This reassures the predictability of our renormalization 
approach, that moreover allows to extend the analysis beyond the ladder approximation 
as we discuss in what follows. 


3.2 Electron self-energy corrections to the ladder approximation 

The most relevant way to improve the ladder approximation corresponds to including 
electron self-energy corrections in the diagrams encoded in the equation of Fig. 2. With 
these additional contributions we may account for an important feature of the electron 
system, which is the growth of the Fermi velocity in the low-energy limit. This leads to a 
reduction of the effective interaction strength, from which we can expect the need to push 
the nominal coupling to larger values in order to reach the phase with dynamical mass 
generation. 

A fully consistent approach can be devised by considering the self-energy contributions 
in the same ladder approximation discussed in Sec. II. In this case we know that the electron 
self-energy coincides with the lowest-order result given by Eq. (2.9). This can be translated 
into a redehnition oi vp-, leading in the fermion propagator to an effective Fermi velocity 

vf(^) =vf + (3-37) 


with 


B{e) 


1 r(l6)F(2-f)F(l-f) 

(47r)i-e/2 F (f - e) 


(3.38) 
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The electron self-energy corrections are then accounted for automatically after replacing 
the parameter vp by in the self-consistent equation (3.15) for the mass vertex. 

It can be easily checked that the perturbative expansion of l/vpi}^ introduced in that 
equation corresponds to the iteration of self-energy rainbow diagrams correcting the fermion 
propagators in the original ladder approximation. 

As in the previous subsection, we can assume that the mass vertex admits now the 
expansion 

CX) 

r™(0,0;k,a;fc) = l +(3.39) 

n=l ' ' 

We can apply a recurrent procedure as before to compute the different orders in (3.39), 
with a main difference in that now each tn is going to depend on all the lower orders in 
the expansion. This is so as the n-th order, when inserted at the improved right-hand-side 
of (3.15), can give rise to contributions to any higher order as the factor l/vF(\i) is also 
expanded. At the end, we arrive at the result that 

n 

tn+i = An+i{e) ^ (-.B(e))”“' ti (3.40) 

1=0 

We have to check again that the vertex can be made finite in the limit e —>• 0 by a 
suitable multiplicative renormalization. In this case, we have first to express all quantities 
in terms of the renormalized Fermi velocity vr arising after subtraction of the pole at the 
right-hand-side of (3.37). vr is defined by 

VF = ZyVR (3.41) 


where has already appeared in (2.11). The renormalized coupling is given now by 


A = 


(3.42) 


4fvr 

When expressed in terms of A, the new renormalization factor Z^ must have the structure 

d*(A) 


Zm — 1 + ■ 


(3.43) 


i=l 


The first terms in the expansions of the residues di{X) can be obtained analytically. 
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with the result that 


di{X) 


d2{X) 

MX) 

d4(A) 


1 1 ,2 76 - 57r2 1908 - 630C(3) - 9l7r2 

"" - ^ - 

168(386 - 91C(3)) - 40047r2 - 337r^ 5 
583207r5 

357r2(615C(3) - 3014) - 15(28028C(3) + 9765C(5) - 94464) + 6327r^ 

5248807r6 

J_ ^ a3 ^ 4 , 14876 - 4410C(3) - 7377r^ 5 

67r2 8l7r3 6487r^ 583207r^ 

604904 - 141204C(3) - 385627r2 - lOOvr"^ 

+- 

1 ^3 1 ^4 218 - ISvr^ ^5 34492 - 11970C(3) - 16297r2 

547r3 1627r4 145807r5 10497607r6 

1 4 , ^ .5 , 1792 - 1357r^ g 

2167r4 8107r5 5248807r6 

^ \5 I ^ \6 I 

6487r5 ^ 32407r6 ^ ^ ' 

A® + ... 

116647r® 


■ A® + .(.3.44) 

(3.45) 
A®+ ..(.3.46) 

(3.47) 

(3.48) 

(3.49) 


The important point is again that the functions di{X) can be chosen with no nonlocal 
dependence (in fact, with no dependence) on the external momentum of the vertex, which 
means that the renormalization can be accomplished with the redefinition of a finite number 
of purely local operators. 

The expansion of the residue di(A) allows us to estimate the effect of the self-energy 
corrections on the anomalous dimension of the vertex. This is given as before by Eq. (3.30), 
while now we have to account for the change in the scaling of the renormalized coupling A. 
This can be obtained from Eqs. (2.13) and (2.14), with the result that 

fi ——A = —eA + ——A^ (3.50) 

a/i 37r 

The computation of the anomalous dimension proceeds then as 


fi dX dZm 
Zm dfj, dX 


(3.51) 


Eq. (3.51) provides an expression for 7 ^ that contains poles of all orders in the e 
parameter. As in the previous subsection, it is therefore highly nontrivial that a finite 
result may be obtained for the anomalous dimension in the limit e —?■ 0. Erom (3.51), the 
equation that we have to inspect in this case is 




\ 2 = 1 / 2=0 2=1 

The zeroth order in e leads to the equation for the anomalous dimension 

7m = 


(3.52) 


(3.53) 
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Figure 4. Plot of the absolute value of the coefhcients in the expansion of di(A) as a power 
series of the renormalized coupling A. 


However, this derivation makes sense only when the cancellation of all the poles in (3.52) 
can be guaranteed, which implies the set of equations 

^Ji+i(A) - Ji(A) ^Ji(A) - I^A (A) = 0 (3.54) 

It is reassuring to see that the analytic expressions given in (3.44)-(3.49) satisfy the 
hierarchy of equations (3.54). A more comprehensive check can be performed however with 
the numerical computation of the expansions of the first residues, that we have carried out 
up to order A^^ and with a precision of 40 digits. In this way, we have been able to certify 
that the conditions (3.54) are also fulfilled at that level of approximation. 

The residue di{X) has in particular an expansion 

OO 

di(A) = ^dS”^A” (3.55) 

n=l 


with coefficients that have been represented in Fig. 4 up to order n = 30. We observe that 

"(n) 

the d\ ' series approaches a geometric sequence at large n. The ratio between consecutive 
orders can be fitted indeed with great accuracy by a behavior like (3.35), allowing to 
compute the asymptotic value 

^ d (3.56) 

in the limit n ^ oo. We have obtained in this way the finite radius of convergence of the 
perturbative series 

Xc = -~^ 1.8660 (3.57) 

d 

with an error estimated to be (as in (3.36)) in the last digit. 

As remarked before, the critical coupling Ac corresponds to the point at which the 
anomalous dimension 7 ^ diverges. This is in turn the signature of the development of a 
nonvanishing expectation value of the mass operator (3.3). Thus we see that, even after 
taking into account the effect of the Fermi velocity renormalization, there still remains a 
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strong-coupling phase in the electron system characterized by the dynamical generation 
of mass for the Dirac fermions. The value of the critical coupling (3.57) is sensibly larger 
than that obtained in the previous subsection, which is consistent with the fact that the 
self-energy corrections effectively reduce the interaction strength. The situation is in this 
respect rather similar to the case of the 2D Dirac semimetals, where the effective growth of 
the Fermi velocity at low energies has been invoked as the reason why no gap is observed 
in the electronic spectrum of graphene[33, 34], even in the free-standing samples prepared 
at very low doping levels about the Dirac point [17]. 

4 Large-approximation for 3D Dirac and Weyl semimetals 

We deal next with an approach complementary to the ladder approximation, paying at¬ 
tention to the effect of the photon self-energy corrections on the electron quasiparticles. 
This will take into account the renormalization of the electron charge, which is a relevant 
feature in the field theory of 3D semimetals as well as in the fully relativistic 3D QED[35]. 
In order to include a consistent set of quantum corrections, we will dress the interaction 
with the sum of bubble diagrams obtained by iteration of the electron-hole polarization, as 
represented in Fig. 5. This approximation can be then considered as the result of taking 
the leading order in a 1/A^ expansion, providing a resolution of the theory in a well-defined 
limit as —>• oo. 



Figure 5. Diagrammatic equations in the large-7V approximation. The thick(thin) straight line 
stands for the dressed(free) Dirac fermion propagator and the thick(thin) wiggly line stands for the 
dressed(free) interaction propagator. 

The corrections to the bare interaction are represented by the polarization n(q,cug), 
from which the full interaction propagator D(q, cug) is obtained according to 

D(q, ujq)~^ = Do(q)"^ - n(q, ujq) (4.1) 

In the present approximation, the polarization is given by the dominant contribution in 
the large-A^ limit 

iIi{(i,ojq) = Nel J -^^^Tr[joGo{p + q_,ujp + ujq)-foGo{p,Up)] (4.2) 
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with the free Dirac propagator 


Go{p,ujp) 


-■youjp + VFl' • p 

—a;2 + Vpp'^ — ir] 


(4.3) 


Computing the integrals in analytic continuation to D = 3 — e, we get the result 

2 2 

u,g) = -NCie) ^^2 _ ^2/^2 )./2 

with 

C(6) = (47r)-/^r (^6) (4.5) 

The polarization (4.4) diverges in the limit e —)• 0, which points at the need to renor¬ 
malize the scalar field 4> mediating the Coulomb interaction. As already mentioned, the 
gauge invariance of the action (2.6) implies that = 1, which means that the electron 

charge is consequently renormalized. In the present approximation, these effects can be 
taken into account at once in the large-expression of the electron self-energy 




iS(k,a;fc) 


d^p dojp 
(27r)^ 2tt 


70 



-70(^fc - ^p) + VFl • (k - p) 

{ujk - uJpY -h u|,(k - p)2 - irP^ 


(l + iVC(e) 


e 


2 

0 


2tz'^vp ( p 2 



(4.6) 


Thus, we can determine the electron charge renormalization by devising a finite limit of 
the effective e-e interaction in (4.6) as e ^ 0. This can be achieved by absorbing the pole 
in (4.4) into the bare charge cq according to the redefinition 


1 _ N 1 

Cq Q'k'^vf e 

In terms of the renormalized charge e, the electron self-energy becomes 


(4.7) 


iS(k,Wfc) 


d^p dojp 


JojuJk - ^p) + VFl • (k - p) 
-(wfc - Wp)2 + vj,(k- p)2 - ir] 



p e 



JVe£_l I JVe^ _ 

e 1 ' 2'It‘^Vp ( p 2— 


(4.8) 


The renormalization program has to be completed anyhow by accounting for the di¬ 
vergences that arise when performing the integrals in (4.8). In order to make sense of the 
theory in the large-A^ limit, we can assume formally that Ne^/2 tt'^vf ~ 0{N^). The dif¬ 
ference between vf and its renormalized counterpart vr is then a quantity of order ~ 1/A^, 
and the electron self-energy gets a natural dependence on the effective coupling 


iVe^ 

^ 2f‘^vr 


(4.9) 
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We may actually resort to a perturbative expansion 


S(k,ia;fc) 


2'k‘^vf e f d^p duip i'jojujk - iOp) + vfI • (k - p) J_ 

N ^ J {2f)^ 2f u|,(k - p)2 + (wfc-a;p)2 p2 


^(_l)n^n+l 

71=0 


1 






(4.10) 


It may be seen from (4.10) that, to order g'^, the self-energy can develop divergences as 
large as ~ 1/e”. If the theory is renormalizable, it must be possible to absorb these poles 
into the definition of suitable renormalization factors in the full fermion propagator (2.10). 
We look then for the large-limit of and Zy, which must have the general structure 


Zib — 



Zy — 



(4.11) 

(4.12) 


with coefficients 5* and Ci depending only on the renormalized coupling g. 

It is not difficult to compute the integrals that are needed to get in general the n-th 
order of S(k, fZJfe). We are interested in terms that are linear in and k, and the only 
technical point is that these contributions must be also regularized in the infrared, using 
for instance the own external frequency or momentum, or some other suitable parameter. 
It can be checked that the choice of a particular infrared regulator is not important when 
extracting the high-energy divergences as e —)• 0. For this reason, we have resorted to 
introduce a hctitious Dirac mass z/, which has the ability of keeping well-behaved both 
types of contributions linear in and k. We rely then on the results for the generic 
integrals 


In = VF 


d^p diOp 
1 


i'joiujk - Wp) 


noWfc 


u|,(k - p)2 -H (Wfc - UJpY + p2 (p2 + I^2j^2^^nel2 

r(^e)r(l-^e) 1 


ne 


(47r)2-^/2 1 - 


r(2-f) 


^|(7i+l)e 


(4.13) 


and 


f d^p doJp vf'J • (k — p) 1 1 

J 2 f u^(k - p)2 -K (aJfc - aJp)2 + p2 (p2 + aj2/.y2 )ne/2 

^ l + (l-e)(l + te)r(^e)r(l-^e) 1 

(47r)2-^/2 (i_e)(3-e) T (2 - f) 


With the help of (4.13) and (4.14), one can obtain analytically for instance the first 
terms in the expansions of the residues in (4.11) and (4.12), by imposing that the renor- 
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malized fermion propagator becomes finite in the limit e —)• 0. We find in this way 


ci(ff) 


C2{9) 

csig) 

CA{g) 

C5{g) 


'24^ 


7 


162 


+ 


5184 


1 


+ 


C(3) \ 


6480 6480 / 


TT 


279936 2799360 


+ 


1 3 

- - 

4_ 

432 ^ 

1 5 

1620 ^ 


1 

1 

2430 

1 


1 

3888 


C(3) 

34992 

/- 


g^ + 


1 


C(3) 


23328 23328 


+ 


69984 


/ + 


8748 


/ + 


1 


5832 


/ + 




and for the Fermi velocity renormalization 


(4.15) 

(4.16) 

(4.17) 

(4.18) 

(4.19) 


h{g) 


h{g) 


h{g) 

h{g) 

h{g) 

beig) 


i i 1 1 2 


1 


9720 

5 


+ 


TT 


+ 


1 

’ ^ 
C(3) 


1 ^ C(3) 


1728 1296 


583200 19440 


+ 


vr 


I C(3) ^ C(5) 


279936 8398080 69984 23328 


-C2{g) - 3 ^ 9 ^ 


1 


1 


1296 


/ + 


1 ^ C(3) 


6480 4860 


TT 


34992 


+ —^1 + 
2099520 69984 J ^ 


—Cs(q) - q - 

jvy; y ^296 


C(3) \ 
i ) 
1 


4860 


1 


+ 


C(3) \ 


23328 17496 / 


/ + 


—cAg) - g'^ 

4vy; 324 y 

1 

-CT,{g) - 


1 


1 


4860 


17496 




1 


1215 


17496 




4374 


/ + 


(4.20) 


(4.21) 

(4.22) 

(4.23) 

(4.24) 

(4.25) 


The important point is that these expansions of the residues do not show any dependence 
on the auxiliary scales ^ and v (or on the external frequency and momentum when these are 
used alternatively to regularize the self-energy in the infrared). This is a nice check of the 
renormalizability of the theory which guarantees that, at least in the large-approxima¬ 
tion, the high-energy divergences can be absorbed into a finite number of renormalization 
factors depending only on the renormalized coupling g. 

The renormalization factors and provide important information about the be¬ 
havior of the electron system in the low-energy limit. In the renormalized theory, the Dirac 
fermion field gets an anomalous scaling dimension 7 ^( 5 ), which can be computed as[31] 


g dZ'ijj 


(4.26) 
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This anomalous dimension governs the scaling of the correlators involving the Dirac field. 
When sitting at a fixed-point of the renormalized parameters, we would get for instance 
for the Dirac propagator 

G(sk, suj) « G(k, w) (4.27) 

On the other hand, the renormalized Fermi velocity vr gets also a scaling dependence 
which can be assessed in terms of a function 'jvig) such that 


fJ- dvR 
VR dg 


lv{g) 


(4.28) 


Thus, the knowledge of 7 ^( 5 ) and 'yvig) allows us to inspect the theory in the limit of long 
wavelengths and low energies as /i —)• 0 . 

The anomalous dimension 7 ^( 5 ) can be computed from the dependence of the renor¬ 
malized coupling g on the auxiliary scale g,, taking into account that 


- dg dZ^ 
dg dg 


(4.29) 


The scaling of g can be obtained in turn by differentiating the expression (4.7) with respect 
to g, bearing in mind the independence of eo with respect to that auxiliary parameter. This 
leads to 


d 


g- 


2 2 

e = —ee 


N 


“dg'' 

As already mentioned, the difference between vp and vr can be taken formally as a quantity 
of order ~ 1 /A^, so that we end up in the large-limit with the equation 


d 1 2 

^ TiO 


This expression can be plugged into (4.29) to get 


(4.31) 


77 


1 1 
'^N 



(4.32) 


Working to leading order in the 1/N expansion, we can set = 1 at the right-hand- 
side of (4.32). The term with no poles in that equation leads to the result 

7b = -^5^ci(ff) (4.33) 

This expression of the anomalous dimension makes only sense, however, provided that 
one can certify the cancellation of the rest of terms carrying poles of all orders in the e 
parameter, which implies the hierarchy of equations 

^Ci+i{g) - ^g^Ciig) = 0 (4.34) 

dg 3 dg 

It is remarkable that the power series of the residues given in (4.15)-(4.19) satisfy 
indeed the hierarchy (4.34). Beyond the analytic approach, we have computed numerically 
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Figure 6. Plot of the absolute value of the coefficients in the expansion of ci{g) as a power 
series of the renormalized coupling g. 

the expansion of the first residues Ci{g) up to order with a precision of 40 digits. In 
this way, we have been able to check that the conditions (4.34) are fulfilled at that level 
of approximation, reassuring the consistency of the large-approach for the present field 
theory. 

Another important detail of the numerical calculation of the residues is the evidence 
that the coefficients of each perturbative expansion approach a geometric sequence for large 
orders of the coupling. In the case of the residue ci{g), we have for instance the power 
series 

OO 

= (4.35) 

n=l 

with coefficients that we have represented in Fig. 6 up to the order we have carried out 

in) 

the numerical computation. The behavior observed for the series c) ^ implies that the 
perturbative expansion must have a finite radius of convergence. This can be obtained by 
approaching the ratio between consecutive coefficients according to the dependence 


which provides indeed an excellent fit at large n. We get in this way the estimate 

r « 0.3333333 (4.37) 

where the error lies in the last digit. We find then the singular point for the coupling g 
(the radius of convergence) at 

= i « 3.0 ± 1.0 X 10"'^ (4.38) 

The critical value gc corresponds to a point where the anomalous scaling dimension 
lijiid) diverges. This can be appreciated in Fig. 7, where we have represented that function 
from the results of our numerical calculation. The singularity found at gc has a precise 


) j-H J.III 

Tr=r-|-1-24- 3+'" 

rj ri n 


(4.36) 
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Figure 7. Plot of the anomalous scaling dimension 7 ^(g), multiplied by the number TV of four- 
component Dirac fermions in the electron system. 


physical meaning, since 7 ^ governs the scaling of all the correlators involving the Dirac 
fermion field. Away from a fixed-point in the renormalized parameters, the scaling is not 
so simple as in (4.27), but from that expression we already get the idea that the divergence 
of 7 ^ implies the decay of the Dirac propagator [36]. In the limit g —)• gc-, the singular¬ 
ity in the scaling dimension leads indeed to the suppression of the quasiparticle weight. 
This characterizes a particular form of correlated behavior, which has been identified in 
several other instances of interacting electrons and has led to constitute the class of so- 
called non-Fermi liquids [37-40]. The distinctive feature of this class is the absence of a 
quasiparticle pole in the electron propagator, which leads to a very appealing paradigm 
to explain unconventional properties like those found in the normal state of copper-oxide 
superconductors. 

In order to ensure the relevance of the divergence at i^c, we have anyhow to verify 
that such a singular behavior is not prevented by some other instability in the low-energy 
scaling of the electron system. We pay attention in particular to the renormalization of the 
Fermi velocity, which has a natural tendency to grow in the low-energy limit. The function 
that dictates the scaling in Eq. (4.28) can be obtained from the independence of the 
bare Fermi velocity on the auxiliary scale g: 

{Z,vr) = 0 (4.39) 


The dependence of on g 
get 


Z„ 


g d 
— 

VR dg 


Working to leading order in 
to the result 


comes only from the coupling g so that, using Eq. (4.31), we 


1 A6,+i(5) + -V V --h{g) = 0 


N 


(4.40) 


i=0 


the 1/TV expansion, the term free of poles in Eq. (4.40) leads 


g d 1 d ^ 

VR dg TV dg 


{9) 


(4.41) 
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Figure 8. Plot of the rate of variation 'jvig) of the renormalized Fermi velocity with respect to 
energy, multiplied by the number N of four-component Dirac fermions in the electron system. 


As in the case of the anomalous scaling dimension, one has to make sure however that 
the rest of terms carrying poles of all orders in Eq. (4.40) cancel out identically, which is 
enforced by the conditions 

The expressions of the residues in (4.20)-(4.25) satisfy indeed the set of equations 
(4.42), and a more extensive numerical computation conhrms that this is also the case for 
the power series of the hi{g) evaluated up to order . This analysis also shows that these 
expansions do not lead to any singularity in the range of couplings up to the critical point 
gc- The corresponding function 7 „(fif) obtained from hi{g) is represented in Fig. 8. The 
regular behavior observed in the plot implies that the scaling of the Fermi velocity is not 
an obstacle for the suppression of the fermion quasiparticles as the interaction strength 
becomes sufficiently large to hit the critical point at gc- 

As a last check, we look also for the possible tendency towards chiral symmetry break¬ 
ing and dynamical mass generation in the large-approximation. For this purpose, we may 
analyze the scaling of the vertex for the mass operator represented in Fig. 9. Computing 
in the limit where both frequency and momentum transfer vanish, we get 


rm(0,0;k,a;fc) = l + i 


d^p dujp 
(27r)^ 27r 


70 


7o(a;fc - uip) + vpl • (k - p) 
-(wfc - a;p)2 + vl{k- p)2 - ig 


70 


X- 






(4.43) 


We have to account as before for the renormalization of the charge, which leads to an 
expression in terms of the renormalized coupling g 


h'm(O) 0, k, ito^) — 


1 , 27r2nF e 

1 + 


/ 


d^p dLJp 1 1 

{2-k)^ 2-k n|,(k - p) 2(wfc-u;p)2 p^ 


E(-o"9”+' 


n=0 




(4.44) 
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Figure 9. Corrections to the vertex for the mass operator in the large-TV approximation. The cross 
represents the operator tjj'ip s-nd the thick wiggly line stands for the dressed interaction propagator 
as defined in Fig. 5. 


As already done in the case of the electron self-energy, we can compute the high-energy 
divergences of the vertex in the limit of vanishing k and regularizing the integrals in 
the infrared with an auxiliary mass v in the Dirac propagators. In this way, the different 
terms in the expansion (4.44) can be obtained using the general result 


f d^p dUp 111 

^ J (27r)^ 27r Vpp"^ + (p^ -|- 

1 2 r(^e)r(l-^e) 1 

(47r)2-e/2 1 _ e r (1 - f) 


(4.45) 


We see that the vertex can develop in general divergences of order ~ 1/e"' at the 
level S'" in the perturbative expansion. These have to be reabsorbed into a multiplicative 
renormalization of the type shown in (3.9). The point to bear in mind is that now 
is made of two different factors, coming independently from the renormalization of the 
composite mass operator (3.3) and the Dirac fermion fields in the vertex[31]: 




(4.46) 


The renormalization factor Z ^2 for the mass operator can have the general structure 


Z^2 — 




(4.47) 


Under the assumption that the theory is renormalizable, it must be possible to choose 
appropriately the residues di{g) to end up with a renormalized vertex, finite in the limit 
e — 0 , given by 

rm,ren — Z^Z^2 T^ (4.48) 

We have checked that the vertex (4.48) can be made free of poles, at least up to order 
we have been able to compute it numerically, with a set of residues di{g) that depend 
only on the renormalized coupling g. We have also seen that these functions have a regular 
behavior in the range up to the critical coupling gc where the anomalous dimension 7 .^ 
diverges. This makes clear that the dominant instability in the large-approximation is 
indeed characterized by the suppression of the electron quasiparticles. 
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Of all the residues di{g), the first of them conveys the most relevant piece of informa¬ 
tion, since it is related to the anomalous dimension of the composite mass operator, defined 
by 

_ g dZ^2 
~ z^2 dfi 

Paralleling the above derivation of similar scaling dimensions, one arrives at the result 

The plot of this function obtained from our numerical computation of di (g) is represented 
in Fig. 10. The regular behavior observed in the figure is the evidence that no singularity 
can be expected in the correlators of the mass operator pm, whose magnitude has got to 
be bound to the scaling dictated by the anomalous dimension 7.^2 in the low-energy limit. 




Figure 10. Plot of the anomalous scaling dimension 7 ^ 2 ( 3 ), multiplied by the number N of four- 
component Dirac fermions in the electron system. 


5 Conclusions 

We have studied the development of the strong-coupling phases in the QED of 3D Dirac 
and Weyl semimetals by means of two different nonperturbative approaches, consisting in 
the sum of ladder diagrams on one hand, and taking the limit of large number of fermion 
flavors on the other hand. We have benefited from the renormalizability that the theory 
shows in both cases, which makes possible to render all the renormalized quantities in¬ 
dependent of the high-energy cutoff. Thus we have been able to compute the anomalous 
scaling dimensions of a number of operators exclusively in terms of the renormalized cou¬ 
pling constant, which has allowed us to determine the precise location of the singularities 
signaling the onset of the strong-coupling phases. 

We have seen that, in the ladder approximation, the 3D Dirac semimetals have a critical 
point marking the boundary of a phase with chiral symmetry breaking and dynamical 
generation of mass, in analogy with the same strong-coupling phenomenon in conventional 
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QED[20-27]. We have found however that such a breakdown of symmetry does not persist 
in the large-A^ limit of the theory, which is instead characterized by the growth of the 
anomalous dimension of the electron quasiparticles at large interaction strength. The 
picture that emerges by combining the results from the two nonperturbative approaches 
is that chiral symmetry breaking must govern the strong-coupling physics of the 3D Dirac 
semimetals for sufficiently small N, while there has to be a transition to a different phase 
with strong suppression of electron quasiparticles prevailing above a certain value of N. 

With the results obtained for the different critical points we can draw an approximate 
phase diagram in terms of the number N of Dirac fermions and the renormalized coupling 
g dehned in (4.9). In principle, the critical coupling (4.38) can provide a reliable estimate 
for the onset of non-Fermi liquid behavior at sufficiently large N. On the other hand, the 
critical value of g deriving from (3.57) may lead to a sensible map of the phase with chiral 
symmetry breaking as long as its magnitude does not become larger than that of (4.38). 
The resulting phase diagram of the 3D Dirac semimetals, represented in Fig. 11, shows the 
regions where the behavior of the phase boundaries may be captured by our alternative 
approximations, away from the intermediate regime about = 4 where the competition 
between the two strong-coupling phases cannot be reliably described within our analytic 
framework. 



Figure 11. Phase diagram of the QED of 3D Dirac semimetals showing an approximate map of the 
phases corresponding to chiral symmetry breaking (CSB) and non-Fermi liquid behavior (NFL), 
obtained from the values of the respective critical couplings in the ladder approach and the large-N 
approximation. 

It is interesting to compare at this point the phase diagram in Fig. 11 with that 
obtained with the resolution of the Schwinger-Dyson equations, which can be trusted for 
all values of N. We can see from the results of Ref. [34] that such a numerical approach 
sets indeed the interplay between the phases with chiral symmetry breaking and non-Fermi 
liquid behavior at = 4. For A^ > 4, we observe that the critical line for that latter phase 
is not straight, although it approaches a constant asymptotic limit at large N. The self- 
consistent resolution leads also to a critical line for chiral symmetry breaking in the regime 
N < 4: that has an approximate linear behavior as a function of A^. It is worth to mention 
that the values of the critical couplings found in the numerical resolution of Ref. [34] are 
at hrst sight much larger than their counterparts in the diagram of Fig. 11. We have 
to bear in mind, however, that the critical values in that paper were given for the bare 
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couplings, in the theory with a finite high-energy cutoff, while critical points like (3.57) 
and (4.38) are referred in the present context to the renormalized couplings. Overall, we 
may conclude that there is good qualitative agreement between the location of the phases 
in the diagram of Fig. 11 and in the more complete map obtained from the resolution of 
the Schwinger-Dyson equations. 

We end up remarking that our analysis can be applied to characterize not only the 
strong-coupling regime of the 3D Dirac semimetals, but also of the Weyl semimetals. In 
these materials, each Weyl point hosts a fermion with a definite chirality, which means 
that the electron system cannot undergo chiral symmetry breaking through condensation 
of some order parameter at zero momentum[41]. This obstruction does not hold however 
for the strong-coupling phase that we have identified in terms of the suppression of fermion 
quasiparticles. It is clear that the large-approach of Sec. IV applies equally well for 
Weyl and Dirac fermions, so that Weyl semimetals are susceptible of developing the phase 
with non-Fermi liquid behavior that we have mapped at large N in the diagram of Fig. 11. 
In this respect, there should be good prospects to observe such an unconventional behavior 
in present candidates for Weyl semimetals, like TaAs or the pyrochlore iridates, which have 
up to 12 pairs of Weyl points[5, 6]. These considerations show that the strong-coupling 
phases studied here are not beyond reach, and that they may be actually realized in 3D 
Dirac or Weyl semimetals with suitably small values of the Fermi velocity or with the large 
number of fermion flavors already exhibited by known materials. 
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